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£N| ■ Abstract 

A classification of gravitating Yang-Mills systems in all dimensions is presented. These systems are set 
up so that they support finite energy solutions. Both regular and black hole solutions are considered, the 
former being the limit of the latter for vanishing event horizon radius. Special attention is paid to systems 
necessarily involving higher order Yang-Mills curvature terms, along with the option of incorporating 
higher order terms in the Riemann curvature. The scope here is restricted to Einstein systems, with or 
without cosmological constant, and the Yang-Mills(-Higgs) systems. 

a -1 

bJO! 1 Introduction 

By gravitating Yang-Mills (YM) fields we understand YM fields on curved backgrounds whose dynamics 
■ includes the backreaction of the gravitational field on it. These are the particle-like and black hole solutions 

fS| \ to the Einstein- Yang-Mills (EYM) systems. 

Initially, EYM solutions [TJ [2] in 3 + 1 dimensional space with Minkowskian signature were primarily of 
interest because they presented configurations with non-Abelian Ziaifl More recently however, gravitating 
non-Abelian solutions have found extensive application in string inspired theories, e.g., in various supergrav- 
ity and D— brane models. These results point to the physical relevance of classical gravitating non-Abelian 
solutions. Moreover in this context, it is mainly such classical solutions in dimensions higher that 3 + 1 that 
find applications. It is our aim here to review gravitating non-Abelian YM solutions in higher dimensions^, 
together with their four dimensional counterparts. 

Soon after the original particle-like (regular) solutions of pQ were constructed, the corresponding black 
hole solutions were found [H [6]. Not long after the discovery of the EYM solutions, both regular and 
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black hole gravitating monopole solutions of the EYM-Higgs (EYMH) systems were constructed in [7] [HI IS] • 
EYM solutions in the presence of a cosmological constant A were constructed later, with negative A in 



pTOl [TT1 [12] , and for A > in [13] . These results indicate that the solutions of the Einstein equations coupled 
to non-Abelian matter fields possess a much richer structure than in the better known U(l) case. They also 
show that our intuition based on solutions with linear field sources may fail in more general situations. 

In this review, our description of gravitating YM systems will in addition to the EYM fields also include 
the Higgs fields. The study of EYM-Higgs (EYMH) systems enables a more extensive description of physical 
phenomena as a result of the symmetry breaking mechanism which (gauged) Higgs models describe. This is a 
result of the special nature of Higgs fields here as dimensional descendents of gauge fields, as will be explained 
below. Therefore we shall restrict our considerations to such Higgs multiplets that afford topological stability 
to the solutions. 



1 These were preceded historically by gravitating Skyrmions 3 which exhibited Skyrmion hair, but the latter have not been 
studied as intensively since. 

2 It is concievable that gravitating Skyrmions may also be of relevance to field theories in higher dimensions, but to date 
there has been no work on them reported in the literature. 
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EYM solutions in (higher) d > 4 dimensions were considered relatively recently. Here two main pos- 
sibilities have been included so far in the literature, which are distinguished by the different asymptotic 
structure of the spacetime. In the first case, of main interest here, the metric is asymptotically Minkowski 
(or (anti-)de Sitter, if a cosmological constant is considered in the action). Such solutions provide the natural 
counterparts of the d = 3 + 1 non-Abelian configurations mentioned above. In the second case, a number of p 
codimensions are included. Such configurations are important if one posits the existence of extra-dimensions 
in the universe, which are likely to be compact and described by a Kaluza-Klein (KK) theory. All known 
KK non-Abelian solutions have no dependence on the extra-dimensions, i.e. these are frozen, as in the case 
of the z— coordinate of the Abrikosov-Nielsen-Olesen vortex. 

A number of results in the literature show that, in the absence of codimensions, the mass of gravitating 
non-Abelian solutions which asymptote to a Minkowski spacetime, is infinite [T3] (the same results holds 
for asymptotically (anti-)de Sitter solutions 1 5j ) . This is not surprising because the usual EYM system in 
d > 5 dimensions does not have the requisite scaling properties for there to exist finite energy solutions. 

However, the scaling properties of the usual EYM system can be altered by the addition of higher 
order terms of the Yang-Mills curvatures. Such terms can occur in the low energy effective action of 
string theory [16], [17]. The hierarchies of both YM and of gravitational (i.e., Einstein) systems will be 
defined below. Employing suitably defined EYM systems featuring higher order YM curvature El terms such 
that the equations of motion remain second order, finite mass/energy static spherically symmetric non- 
Abelian solutions in d > 4 dimensions were constructed. The existence of this type of configurations is a 
nonpcrturbative effect, since they cannot be predicted in a perturbative approach around the solutions of 
the usual EYM system. 

With zero cosmological canstant, d > 4 regular solutions of the extended EYM system were found 
and analysed in jT5J QI5] , and black hole solutions in [2U] , [2T] . Both regular and black hole with negative 
cosmological canstant A were presented in [2T]. With positive A likewise, such EYM solutions were given 
in [25]. Thus there is a comprehensive sample of finite mass, static, spherically symmetric EYM field 
configurations in d = D + 1 dimensional spacetimes. These are the higher dimensional counterparts of the 
EYM solutions in four dimensions, and in general their non-Abelian matter content is composed of several 
YM terms with various (appropriate) scaling properties. Such models proliferate with increasing dimension 
and exhibit additional features, absent in the usual 3 + 1 dimensional EYM casfl 

Concerning EYMH solutions describing gravitating monopoles in higher dimensions, recently a very 
particular hierarchy of YMH models in d = Ap dimensions has been studied where regular and black hole 
solutions are constructed [24]. The reason for restricting to a particular family of models is the ubiquity of 
Higgs models in dimensions higher than 3 + 1, as will be described below. 

All the above noted EYM and EYMH solutions are static and spherically symmetric. But in 3 + 1 
dimensions there are axially symmetric solutions [25 ,26,27 with very interesting properties, so it is natural 
to seek their higher dimensional counterparts. To the best of our knowledge, the only result in this direction 
is that in [28 , for the simplest EYM system in d = 4+ 1 dimensions. The symmetry imposed in that case 
was bi-azimuthal symmetry in the 4 spacelike dimensions. 

Concerning the case of non-Abelian solutions with codimensions, the situation is less explored, the case 
d = 5 being the only one discussed in a systematic way in the literature. Both non-Abelian vortices and 
black strings with non-Abelian hair have been studied in the literature, starting with the pioneering work of 

M- 

Finally, we note that all the work referred to above pertains to gravitating YM in D + 1 dimensional 
spacetime with Minkowskian signature. This is because to date all of the work on Euclidean EYM in all 
dimensions is effectively restricted to the study of YM fields on fixed (gravitational) backgrounds. This, even 
though the earliest EYM solutions were constructed in 4 Euclidean dimensions (see [2^1 [30] an d [HUGH]). As 
such, these results are not central to the present review which is mainly concerned with fully backreacting 

3 Higher order Riemann curvature terms cannot be employed for this purpose for when they are of sufficiently higher order 
to satisfy the scaling requirement, they invariably vanish or are total divergence terms. 

4 The exception is a very particular hierarchy of gravitating YM models in d = 4p dimensions, consisting of a single higher 
order YM curvature term which scales appropriately. In this case |23j all the qualitative features of the EYM solutions in d = 4 
are preserved. 
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matter-gravity solutions. 

The outline of this paper is as follows. In the next section, 2, we introduce the hierarchies of gravitational 
and Yang-Mills (and Higgs) systems, in all dimensions. Then in the following section, 3, we review known 
results in gravitating YM and YMH solutions in higher dimensions, including a brief description of EYM 
solutions in Euclidean space. There we pay special attention to the existence and the generic properties 
of the non-Abelian solutions and not to their physical applications. In section 4 we give a summary and 
outlook. 

2 Einstein and Yang— Mills hierarchies 

The gravitating solutions of interest are those with finite energy in the Minkowskian regime, and with finite 
action in the Euclidean case. They would describe spherically symmetric black hole solutions of horizon 
radius ru, which include the regular particle like configuration in the limit of ru — 0. As noted in section 
1, the non-Abelian matter sector will consist both of YM, and YM-Higgs (YMH) fields. These systems are 
defined on Euclidean spaces since in the Minkowskian case the solutions in question are static so that the 
fields depend on the spacelike coordinates only. 

Even before gravitating the YMH field, the existence of finite energy /action solutions is contingent on the 
requisite scaling requirement being satisfied. In the case of flat (Euclidean) space, this amounts to satisfying 
the familiar virial relation, as long as the (static) Hamiltonian/Lagrangian is positive definite. In the case 
of a gravitating system the property of positive definiteness is absent, since all Einstein systems are defined 
such that they are not bounded from below. One can procede nonetheless heuristically seeking to satisfy 
the virial relation. A systematic analysis of this for higher dimensional EYM systems was presented in |21j . 
and we will return to this question in section 3 when the static systems are gravitated and will elaborate 
on the scaling arguments. For now, we restrict our attention to the YMH systems on a flat backround. 
It is the virial constraint that necessitates the introduction of higher order curvature terms in the (static) 
Hamiltonian/Lagrangian. These are the members of the Einstein and YM hierarchies to be defined below. 

There is a marked difference in the status of the Einstein hierarchy and the YM hierarchy in the present 
context. The higher order YM terms are necessary for rendering the scaling properties of the EYM system in 
question appropriate for supporting finite energy solutions, while the correspoding terms in the gravitational 
hierarchy do not play this role. All the higher order gravitational terms {e.g., Gauss-Bonnet) with the 
required scaling property either vanish or are total divergence. 

Anticipating the definition of the Einstein hierarchy in the next subsection, we note that the p— Einstein 
system eRu,) as defined by the relations |T]), or ([3|) bellow, scales as L~ 2p . (Note that the usual Einstein- 
Hilbert term,( or 1— Einstein system in the terminology of this work), scales as L~ 2 and the usual YM, or 
1— YM system, scales as L~ A .) According to the heuristic scaling argument, adding the p— Einstein term to 
the usual EYM system would result in the correct scaling if d < 2p + 1 . In the limiting case when d = 2p, 
e R( v ) is a total divergence and beyond that it vansihes, as will be seen from the definitions below. 

On the other, hand higher order members of the gravitational hierarchy play a quantitatively interesting 
role in highlighting certain qualitative features of the solutions, that repeat in dimensions modulo 4p. They 
are also of intrinsic interest in some considerations for the case of Euclidean signature. For those reasons, 
they are included in this review. 

In the next subsection, 2.1, we present the definition of gravitational systems in all dimensions, which 
in this review we refer to as the hierarchy of Einstein systems. In general, higher dimensional gravitational 
systems are composed of the superposition of individual members of the Einstein hierarchy. In subsection 
2.2 we define the Yang-Mills (YM) systems in all even (spacelike) dimensions, referred to as the hierarchy 
of YM systems here. YM systems in all dimensions, both odd and even, can then be constructed from the 
superposition of individual members of the YM hierarchy. Then in subsection 2.3 we introduce that subclass 
of (gauged) Higgs models that are relevant to the presentation here, namely those YM-Higgs (YMH) models 
that have been gravitated to date, along with an example of the next most natural candidate. Finally in 
subsection 2.4 we state the expressions for the reduced Lagrangians of the systems introduced, subject to 
the appropriate symmetries. This will be mainly the case of most interest involving the static fields which 
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are spherically symmetric in the spacelike dimensions. 



2.1 Gravitational systems in all dimensions: Einstein hierarchy 

The gravitational systems which we shall refer to as the Einstein hierarchy in d = 2p+k spacetime dimensions 
are defined as the product of the determinant of the Vielbeine e™, e = det e™, with Ricci scalars R( p ), defined 
by 

r ^ — p R, s — <?I*1H2—H2pviva...i'k f,ni n 2q n k r>m 1 m 2 ...m 2p r-i\ 

*-fi-E — — t e 1/1 C V2 ...c Uk t-m 1 m 2 ...m 2p n 1 n 2 ...n k lx -fi lP , 2 ...M2p ' V / 

R^i^. '/lTp 2 '' = R(2p) being the p-fold totally antisymmetrised product of the Riemann curvature, in compo- 
nent notation 

R(2) = RTu 12 = d lu , to" 1 ]" 12 + lu™ 1 " tu n ™ 2 , 

v ' M1M2 [Mi M2J [Ml M2J ' 

w™' 1 being the Levi-Civita spin connection. It is clear from the definition |T|) that the spacetime dimen- 
sionality d = 2p + k sets an upper limit on the highest order nontrivial member of this hierarchy, the term 
R(p)(k=o) being the (total divergence) Euler-Hirzebruch density. We shall refer to ei?( p ) in (Q} asp— Einstein 
systems (i.e. p = 1 is the usual Einstein-Hilbert Lagrangean, p = 2 is the Gauss-Bonnet term etc.). 

Subjecting fl} to the variational principle with respect to the arbitrary variation of the Vielbeine one 
arrives at what we refer to as the p-Einstein equation 

G (P)M = R (p)7 ~ R (P) e M ' ( 2 ) 

in terms of the p-Einstein tensor G( p )™, with R^ and R( p )™ being the p-th order Ricci scalar and the p-th 
order Ricci tensor defined respectively by 

R, ■. pm 1 m 2 ...m 2p Ml P Ml P M2 P (o\ 

1 HP) Jl Ml H2---U2n c mi c m 2 ■ ■ ■ c m 2p \° ) 



D m r>mrn 2 m 3 ...m 2v n 2 p, 3 fi 2p (a\ 

Jt (p)M L MM2M3---M2 P c m 2 c m 3 • • • c m 2p ' 



The p— Einstein hierarchy of gravitational systems is defined by e times the p— Ricci scalar (|3]), and the 
most general Einstein system is given by the maximal number of nonvanishing superpositions of all p— Ricci 
scalars. 

An interesting property of the p— Riemann and q— Riemann curvatures in even dimensions is the doublc- 
sclfduality constraint 

pmim 2 ..m 2p l_ ,„2(p-q) e T>v^v 2 ....v 2q m 1 m 2 ...m 2p nin 2 ...n 2q fr\ 

J VlM2----M2p — [(2q)!] 2 f * 1 ^---V2 p VlV 2 ...V 2 , ] il mtl 2 ..J»2, fc > \°J 

where the Hodge dual of the 2q— form curvature is equated to the 2p— form curvature, with the dimensionful 
constant k compensating for the difference in the respective dimensions. ([5]) can be stated for both Eu- 
clidean and Minkowskian signatures, with the ± sign respectively.For Euclidean signature, contracting the 
constraint ([1J with the appropriate number of vielbeine one arrives at the vacuum Einstein equations for the 
(p, q)— Einstein system 

£ (p , g) _ E -e(i? (p) + K 2 ^i? (9) +A), (6) 

where Ri v ) and Rr q \ in ([6]) are defined by ([3]) and A is a cosmological constant whose value is related to the 
constant k, as shown in [33]. This is not valid in the case of Minkowskian signature (see Appendix of [34] 
for details). 



2.2 The Yang— Mills hierarchy 

Since we seek finite energy/action solutions of the gravitating YM systems, the relevant members of the YM 
hierarchy to be defined, are those which support finite action solutions in the spacelike (Euclidean) subspace 
D of the spacetime d = D + 1. The flat space solutions in question will have topological stability when D is 
even. To avail of topological stability for these finite energy/action solutions in odd dimensions D, one has 
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to consider the Higgs models derived from the dimensional descent of the YM systems in even D + N with 
compact odd codimension N, The YM hierarchy is presented in the section 2.2.1 and the ensuing YMH 
models in section 2.2.2, below. 

Using the notation F{2) — for the 2-form Yang-Mills (YM) curvature, the 2p-form YM tensor 

F{2p) — F{2) A F(2) A ... A F{2) , p - times (7) 

is a p fold totally antisymmetrised product of the 2-form curvature. 

The p—YM system of the YM hierarchy is defined, on R 4p , by the Lagrangian density 

TrF(2p) 2 . (8) 

In 2n dimensions, partitioning n as n = p + q, the Hodge dual of the 2g-form field F(2q), namely 
(*F(2q))(2p), is a 2p-form. 
Starting from the inequality 

Tr[F{2p) - k *F{2q)] 2 > , (9) 

it follows that 

Tr[F(2p) 2 + k 2 F(2q) 2 } > 2k C n , (10) 

where C n is the n-th Chern-Pontryagin density. In ([9]) and (|10p , the constant k has the dimension of length 
to the power of (p — q). 

The element of the YM systems labeled by (p, q) in (even) 2(p + q) dimensions are defined by Lagrangians 
defined by the densities on the left hand side of (|10p . When in particular p = q, then these systems are 
conformally invariant and we refer to them as the p—YM members of the YM hierarchy. 

The inequality (|10[) presents a topological lower bound which guarantees that finite action solutions to 
the Euler-Lagrange equations exist. Of particular interest are solutions to first order self-duality equations 
which solve the second order Euler-Lagrange equations, when (fit?)) can be saturated. 

For M 2 ™ = R 2 ™, the self-duality equations support nontrivial solutions only if q = p, 

F{2p) = *F(2p) . (11) 

For p = 1, i.e. in four Euclidean dimensions, is the usual YM selfduality equation supporting instanton 
solutions. Of these, the spherically symmetric [351 1M] and axially symmetric [37J [Ml (SHj instantons on R 4p 
are the known. For p > 2, i.e. in dimensions eight and higher, only sphericaly symmetric [36] and axially 
symmetric 38, 39 solutions can be constructed, because in these dimensions (TTTT) are overdetermined [ID] . 

In the r>l region, all these 'instanton' fields on R 2 ™, whether self-dual or not, asymptotically behave 
as pure-gauge 

A -> gdg- 1 

For M 2n = G/H, namely on compact coset spaces, the self-duality equations support nontrivial solutions 
for all p and q, 

F(2p) - k *F{2q) (12) 

where the constant k is some power of the 'radius' of the (compact) space. The simplest examples are 
M 2n = S 2n , the 2n-spheres 01], and M 2n = CP™, the complex projective spaces 02] . 

The above definitions of the YM systems can be formally extended to all dimensions, including all odd 
dimensions. The only difference this makes is that all topological lower bounds enabling the construction of 
instantons are then lost, but this is immaterial from the viewpoint in the present review. 

2.3 Higgs models on R D 

Higgs fields have the same dimensions as gauge connections and appear as the extra components of the 
latter under dimensional reduction, when the extra dimension is a compact symmetric space. Dimensional 
reduction of gauge fields over a compact codimension is implemented by the imposition of the symmetry of 
the compact coset space on the coordinates of the codimensions. In this respect, the calculus of dimensional 
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reduction does not differ from that of imposition of symmetries generally, which is the relevant formalism 
used in this, 2.3, and the next subsection 2.3. 

The calculus of imposition of symmetry on gauge fields that has been used in the works being reviewed 
here is that of Schwarz [43, 44, 45 . This formalism was adapted to the dimensional reduction over arbitrary 
codimensions in HZl 132] • 

In general one can employ a linear combination of inequalities (fTU|) . for all p < D/4 and q < D/4. 
Restricting, for simplicity to the 4p dimensional conformal invariant systems in (|10[) . i.e. to p = q = D/4, 
the descent over the compact space K 4p ~ d is described by 

f T{2p) 2 > f C 2p , (13) 

JR d xK 4 p- d JB, d xK 4 p- d 

where T{2p) is the 2p— form curvature of the 1— form connection A on the higher dimensional space H D x 
K 4p ~ D . Imposing the symmetry appropriate to K 4p ~ D on the gauge fields results in the breaking of the 
original gauge group to, say, the residual gauge group g for the fields on R D . Performing then the integration 
over the compact space K 4p ~ D leads to the Lagrangian C[A, </>], of the residual Higgs model on R D . A here 
is the connection taking values in the algebra of g and <f> is the Higgs multiplet whose structure under g 
depends on the detailed choice of K 4p ~ D , implying the following gauge transformations 

A -> gAg- 1 + gd g- 1 

and depending on the choice of K 4p ~ D , 

<\>^> 9<t> 9~ X , or , (j>^ g<j> , etc. 
The inequality (fT0|) leads, after this dimensional descent, to 

/ > / v-n[A,c/>]= [ n[A,<i>], (u) 

where C[A, <fi] = C[F, Dtp, \4>\ 2 , r] 2 ] is the residual Lagrangian in terms of the residual gauge connection A and 
its curvature F, the Higgs fields <j) and its covariant derivative D(f> and the inverse of the compactification 
'radius' r\. The latter is simply the VEV of the Higgs field, seen clearly from the typical form of the 
components of the curvature F on the extra (compact) space K 4p ~ D 

F\ Kip -o ~ {r, 2 - \(t>\ 2 ) ® E => lim|0| 2 = 77 2 (15) 

r — >oo 

where E are, symbolically, spin-matrices/Clebsch-Gordan coefficients. 

It should be noted at this stage that subjecting the selfduality equations (fTTj) to this dimensional de- 
scent results in Bogomol 'nyi equations on R D , which for p > 2 in R D x K 4p ' D (cf. fl3])) turn out to be 
overdetermined [40] with few exceptions. 

There arise a plethora of Higgs models, depending on the mode of dimensional descent, namely on the 
particular choice of the compact codimension K 4p ~ d . We will not dwell on various modes of descent and the 
detailed properties of the descended YMH models here, and will limit our attention to those models that 
have been gravitated to date. 

Perhaps the most interesting, or useful, family of YMH models on M D arrived at via this descent 
mechanism are those in which the residual gauge group is SO(D), and the Higgs field multiplet is an isovector 
of SO{D) B It turns out that only when D = 2 and when D = 4p — 1 in the descent over 
the resulting Bogomol'nyi equations are not overdetermined [40 . The D = 2 case is uninteresting from 
our present perspective since the Abelian Higgs systems in that case live in 2 + 1 spacetime dimensions and 
gravitating them is unproductive. This leaves the family of SO(D) Higgs models that live in D = 4p — 1 
space, or d = 4p spacetime dimensions, for which the flat space Bogomol'nyi equations can be saturated. 

5 Emlpoying Dirac matrix representations for the algebra of SO(D) in terms of i,j = 1,2, . . . ,, the Higgs field takes its 
values in the matrix basis r^xJ+l- with (ry,rj .D+i) representing the algebra of SO(d) = SO(D + 1) 



G 



These are the only Higgs models to date, that are gravitated [21] • (The flat space solutions of this system 
was studied in [3B], which are direct generalisations of the usual BPS monopoles with p = 1.) 

The YM field in the YMH models on R D discussed thus far, is purely magnetic supporting a 'magnetic' 
monopole. But when it comes to YMH models, as stated earlier, the presence of the Higgs field enables 
the support of a dyon in d = D + 1 dimensional spacetime. In the usual [39] sense as the dyon in 3 + 1 
dimensions, the Higgs field partners the newly introduced 'electric' YM potential A . Thus we can describe 
SO(d) dyons0 in d— dimensional spacetime. 

When the dimension of the spacetime d is even, then the chiral representations of the algebra of SO(d) can 
be employed, namely replacing the Dirac representation matrices r p „ = (IYj, T^) with nffi — (Sy , )> 
the precise definition of these matrices to be given explicitly in the next subsection. 

When by contrast the dimension of the spacetime d is odd, then the Dirac represenation matrices are 
the appropriate ones to be used, e.g., in d = 4 + 1 spacetime |50) when D = 4. (No such higher dimensional 
monopoles are gravitated to date.) 

The family of YMH in d — Ap dimensional spacetime that are gravitated result from the simplest mode 
of descent over n D x K 1 , with D = Ap — 1, i.e., over one codimension. The Lagrangian densities can be 
expressed for arbitrary p in flat space as |48j 



-p— YMH 



= Tr 



= Tr 



F(2p) 2 + 2p (F(2p - 2) A D<S>f 



(Fm 1 m 2 ...m 2p ) + 2p (F[ mim2 .m 2p -2 ^«2 P -l] < ^ > ) 



(16) 



in an obvious notation. 



2.4 Static spherically symmetric fields 

Since almost all the work reviewed in this article involves static spherically symmetric only, we will subject 
the above introduced systems to this symmetry only. 

The usual metric Ansatz with spherical symmetry in d — 1 dimensional subspace is 

ds 2 = TN(r)a 2 (r)dT 2 + Nir^dr 2 + r 2 dQ 2 {d _ 2) , (17) 

the T sign pertaining to Lorentzian and Euclidean signatures, and with dtl? d _ 2 \ being the metric on S d ~ 2 . 

Subject to the Ansatz (fTT|) . the reduced one dimensional Lagrangian of the p— Einstein system (TTJ (or 
([3])), in d— dimensional spacetime is calculated. After neglecting the appropriate surface terms, this can be 
expressed compactly as 

L M *v_ } d - 2 ^ * r r d-ap-i(i_jv)rl . ( i8 ) 

(grav) 2 2 P-! (d - 2p - 1)! dr 1 V ' 

We next state the static spherically symmetric Ansatz for the p— YM system (|8]) and the p— YMH system, 

(1 - w(r)\ 

$ = t) h(r) Xj T hd , A Q = u(r) xj Tj^ d , A; = TyXj . (19) 



for odd dimensional spacetime d, and, for even d 

$ = n h(r) Xj , A = u(r) Xj , A t = (j—^j 4? &i , (20) 

where 



6 While both the Higgs field and Aq take their values in the Dirac matrix basis m the in the 3 + 1 dimensional 

case the 'enlarged' algebra SO (4) splits in the two SU(2) subalgebras, whence the dyon and the monopole are both described 
by SU(2) matrices. In all higher dimensions, this is not the case and the full algebra employed is that is SO(D + 1), where 
d = D + 1 is the dimension of the spacetime. 
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are the chiral representations of SO(d). The dimensionful constant r\ in (fl~9"l) - (f2"01) is the Higgs VEV. 

It should be stated here that the choice of gauge group here is made with the purpose of enabling the 
construction of nontrivial finite energy solutions, and that this is the minimal size of gauge group in each 
case. Larger gauge groups, with the appropriate representations containing these, can be employed in case 
of necessity, e.g., when a Chern-Simons term is to be introduced in the Lagarangian. 

Imposing this symmetry, i.e., substituting (|19p or (|20[) into the appropriate Lagrange density, results 
in exactly the same one dimensional reduced Lagrangians in both cases (except for an overall factor of 2). 
This is because the algebraic manipulations involved in both cases are identical. The situation changes if 
Fermions are introduced, but we do not do that here. 

The resulting reduced one dimensional Lagrangian of the p—YM system (|SJ), augmented by the Higgs 
kinetic term in (|16p . in d— spacetime dimensions is 

fw'\ 2 /l-w 2 
(2p)iV - +(d-[2p + l}) ' 



YMH 2- (2p)\ 1 (d- [2p+l})\ r 



=F 



(d-2)! (2p-l) 
(d-2p)\ a 



r I \ r 

2 



2 



(d-2)! (2p-l) 2 
+ {d-2p)\ a V 



([(1 - W *f-\]'? + (d~ 2 P ) { \^-[(1 - w*f-'uf (H) 

([(i - W ^h\r + (d- 2jp)^=%i - wy- ^ (^) 2 



(21) 



This is the most general matter part consisting of the YMH system for this particular family of Higgs models. 

Before proceeding to describe the various types of gravitating YM and YMH solutions in the next section, 
let us make some remarks concerning the particular choices of the various models employed. 

• To recover the formula used for the p—YM system (|8]), one simply replaces ft, = in the third line of 
(f2~Tj) . thus eliminating the Higgs field. 



To recover the matter Lagrangian used for the gravitating monopoles in d = Ap dimensions presented 
in [24], one replaces u — > in the second line of (f2Tj) . and, sets d = Ap since the gravitating monopoles 
there are constructed only in those spacetime dimensions. Of course, it would be possible also to 
gravitate dyons in higher dimensions just as in d = 4 (see for example [51] and references therein), but 
this has not been done to date. 

The choice of model for the monopoles gravitated in [24] was made firstly such that the YM term 
and the Higgs kinetic term have the same dimensions. This is precisely with the criterion that there 
should not be a mismatch of dimensions giving rise to a conical fixed point sigularity, so as not to 
cloud an otherwise more complicated system. This family of models is dimensionally descended from 
the p—YM system and hence the Bogomol'nyi equations can be saturated, which is unimportant since 
the backreaction of gravity prevents this saturation anyway. Secondly, we opted for gravitating the 
p— Higgs models with p— gravity for the more or less aesthetic reasons of that choice in the EYM case 
of [23]. 

In the d > 5 EYM case, various matter systems consisting of the superpositions of p—YM systems 
(JSJ) are gravitated with the 1— Einstein gravity (usual Einstein- Hilbert Lagrangian). This immediately 
introduces a mismatch between the dimensions of the constituent terms in the Lagrangian. It is found 
that one result of this is the absence of the radial excitations (higher node solutions) observed in the 
p = 1, d = 4 Bartnik-McKinnon [T] case. Another result is that in addition to the Reissner-Nordstrom 
fixed point, there arises a new singularity |20] in d = 4+ 1 dimensions. This singularity was found 19 
to repeat in d = 4p + 1 dimensions, modulo Ap. The fixed point corresponding to it was called a 
conical fixed point. 

The electric potential necessarily vanishes for asymptotically flat finite energy solutions of the EYM 
system, i.e. u(r) — if h(r) — 0. The proof here is similar to that found in [52], [53] for d = 3 + 1 
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dimensions. One starts with the equation for the electric potential u(r), which, for a generic model 
with a number of P terms in the YM ierarchy can be rewritten as 



P i / r d-2 \' P r d-2 f 

p=l ^ ' p=l ^ 



+ 



2(d-2p)(2p- 1) 
iV 



IT" 



(here, to simplify the relations, we denote W = (1 — u? 2 ) p 1 ). One can easily see that the r.h.s. of (|22l) 
is a stricly positive quantity. Thus the integral of the l.h.s. should also be positive, 



(where r$ — 0, r# for particle- like and black hole solutions, respectively). However, the regularity of 
the solutions together with finite energy requirements impose that, in the above relation, both the 
contributions at r = r$ and at infinity vanish. As a result, u(r) should vanish for any reasonable 
solution. The same proof generalises for anti-de Sitter solutions, the only exception being the systems 
featuring exclusively the p-th terms of the hierarchy, in d — Ap dimensions. 

• Departing from these generic models, there is a family of models for which the mismatch of the 
dimensionality of the constituent terms is removed. Like in the usual EYM system consisting of the 
1— Einstein and 1— YM systems in d = 3 + 1, this family of models [23] consists exclusively of the 
p— Einstein and p— YM systems in d = 4p. The result is that all qualitative features of the EYM 
solutions of PQ are preserved. Indeed, if p— Einstein is replaced by q— Einstein, q ^ p, the salient 
features persist but are quantitatively somewhat deformed. 

3 Non-Abelian solutions in d— dimensions 

3.1 Solutions with Lorentzian signature 

3.1.1 Einstein Yang-Mills solutions in four dimensions 

The closed form solutions in Chakrabarti et al [31], [32] are probably the first examples of black holes 
with non-Abelian hair, albeit with no backreaction between gravity and the non-Abelian matter. Fully 
sclfgravitating EYM solutions were constructed somewhat after those on fixed backgrounds, by Bartnik 
and McKinnon [1J. These were regular particle- like solutions, and were soon followed by their black hole 
counterparts in [1] [5] [5] . 

Subsequently, a large literature has developed on this subject, extending to systems with a cosmological 
constant, and separately, to systems whose Lagrangian contains also a Higgs field, supporting gravitating 
monopoles. Extensions of the EYM system to include other fields which enter various stringy models have 
been considered as well, in particular for a Gauss-Bonnet quadratic curvature term coupled with a dilaton 
|54j . However, these solutions are beyond the scope of the present review. (A detailed review of the various 
d = 4 gravitating solutions with non-Abelian fields was presented a decade ago in [2j. The case of d = 4 
asymptotically anti-de Sitter (AdS) solutions which was not covered in [2], was the subject of the recent 
review [55].) 

These were all static spherically symmetric solutions, some of whose salient properties will be contrasted 
in their higher dimensional counterparts to be reported in the next subsection. Restricting to solutions with 
a gauge group SU(2), their basic properties are: 

In the EYM case, the asymptotically flat solutions 

• were sphalerons, i.e. that they were unstable [5B], [57] since there was no topological charge to supply 
the energy with a lower bound, and, 




(23) 
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• they present radial excitations characterised by a number k of nodes of the function magnetic gauge 
w(r) in (j2"0"]k 

• the non-Abelian electric potential necessarily vanishes for both globally regular and black hole solutions 
with finite energy [52], [53] . 

When a negative cosmological constant is added to the EYM Lagrangian [TU], [H], the asymptotically AdS 
solutions exhibit new and interesting features, namely that now 

• the asymptotic value of the function w[r) in (|20p is not fixed a priori, which leads to finite mass 
solutions with a nonvanishing non-Abelian magnetic charge, even without a Higgs field; 

• stable solutions have been shown to exist [TO], [58] (this corresponds basically to the case where the 
profile of the function w(r) presents no nodes); 

• black holes with non-Abelian hair and a nonspherical topology of the event horizon have been found 
for A < in [53], [60]; 

• most importantly in this case, it becomes possible to construct finite energy solutions with a nonvan- 
ishing Aq [11] . i.e. non-Abelian dyons; 

• moreover, finite mass solutions with AdS asymptotics exist for any A < 0. 

In the case of a positive cosmological constant [63], [13], [61], [62], by contrast, 

• the EYM solutions with de Sitter (dS) asymptotics exist for sufficiently small values of A only; 

• all solutions have been shown to be unstable, since w(r) necessarily presents nodes (although the 
asymptotic value of the magnetic gauge potential is not fixed a priori); 

• the electric potential Aq necessarily vanishes for all dS solutions. 

In the EYMH (gravitating monopole) case, which differs from the EYM in that a dimensionful constant (the 
Higgs vacuum expectation value) appears in the Lagrangian, the solutions with A < 

• are topologically stable in the YMH sector, stabilised by the monopole charge; 

• they present radial excitations characterised by multinode profiles in the function w(r) in (|20p as in 
the YMH case, and in addition, 

• due to the presence of the dimensionful constant in the Lagrangian, they exhibit a Reissner-Nordstrom 
fixed point, which results in the absence of solutions for a range of the gravitational coupling constant. 

The picture is more complicated for asymptotically dS gravitating monopole solutions. Refs. [64j . [62] 
presented arguments that 

• although the total mass within the cosmological horizon of the monopoles is finite, their mass evaluated 
at timelike infinity generically diverges; 

• no solutions exist in the absence of a Higgs potential. 

The d = 4 asymptotically Minkowski (or AdS) EYM solutions discussed above have axially symmetric 
generalisations. The first work in this direction was [25j . which presented a generalization of the Bartnik- 
McKinnon solutions characterized by a pair of integers [k, n), where n is an integer - the winding number 
and k is the node number of the amplitude w(r). The black hole counterparts of these configurations were 
discussed in [2BJ, which shows that Israel's theorem [65] does not generalise to the non-Abelian case (i.e. a 
static black hole is not necessarily spherically symmetric). These asymptotically flat solutions were extended 
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afterwards in various directions, see e.g. |66j . [67] , [68] . |69j . They present also generalizations with a negative 
cosmological constant, which were discussed in [70], [7T], [75]. The case of axially symmetric non-Abelian 
solutions with dS asymptotics was not considered yet in the literature. 

Interestingly, although non-Abelian generalizations of the Kerr-Newmann black hole were shown to exist 
[73], [71], it turns out that the Bartnik-McKinnon globally regular solutions admit no asymptotically flat 
rotating generalizations |75j . [73j (however, note that they were predicted in a perturbative approach |76]V 
Not completely unexpected, spinning EYM solitons were found to exist for AdS asymptotics [77], [72]. 

Finally, let us remark that both the EYM and EYMH systems present nontrivial solutions with a NUT 
charge [78] , [79] . These solutions approach asymptotically the Taub-NUT spacetime [80] and provide the non- 
Abelian counterparts of the U(l) Brill solution [81]. The nonexistence results in [52], [53] are circumvented 
by these asymptotically locally flat solutions, which necessarily present a nonzero electric part A$ of the 
non-Abelian potential. 




Figure 1: Typical particle-like higher dimensional EYM-7 ? (2) 2 solutions in asymptotically flat and anti-de Sitter 
spacetimes (left) and de Sitter spacetime (right). The function m(r) corresponds to the local mass-energy density. 



3.1.2 Einstein Yang-Mills solutions in higher dimensions 

Gravitating non-Abelian fields in higher dimensions have been considered for the first time in [14j for d = 
5 and a YM model containinig the usual F(2) 2 term only. For spherically symmetric regular solutions 
asymptoting to the Minkowski background, it was found that their energy is infinite. Then in Ref. [15] 
it was proven that the energy of the black hole is also infinite. Moreover, [TS] extended these results to 
asymptotically AdS solutions. When employing the Einstein-Hilbert gravity only, one usually defines 

N(r) = i-I^0 ^ (24) 

[1 r d -3 (d-2)(d-lY [ ' 

the function m(r) being related to the local mass-energy density up to some d— dependent factor. The results 
in [14] . [15] prove that, in five dimensions, m(r) — > logr, as r — * oo. 

As discussed in [21], this is a generic feature of all higher dimensional EYM solutions with a F 2 term 
only (i.e. Lym = T iF^ v F a ^). Although these configurations are still asymptotically Minkowski, their mass 
function generically diverges as r d ~ 5 (or as logr for d = 5). A similar conclusion is reached when considering 
[2T] solutions of a EYM-A model containinig the usual F 2 term only and aproaching asymptotically an AdS 
(or dS) background^- This can most easily be seen by considering the simplest w(r) = solution of the 



7 Asymptotically AdS solutions with diverging mass have been considered by some authors, mainly for a scalar field in the 
bulk (see e.g. [82]). In this case it might be possible to relax the standard asymptotic conditions without loosing the original 
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EYM equations. This corresponds to the gravitating Dirac-Yang monopoles [SSI [57], which are non-Abelian 
configurations (except in d = 4 where one has the Abelian Dirac monopole) . These fields are singular at the 
origin and hence an event horizon should be present. 

The result is a black hole solution, which for d > 5 has a line element 

9 dr 2 , ~ / u 2 M 2Ar 2 \ , 
ds = 2 n + r 2 dfif d 21 - 1 - ^ - -r\ - — — — dt 2 , (25 

l_i^__Mo 2Ar2 ( d -2) I r 2 r d-3 (d - 2)(d - 1) J V ' 

where /i is a constant fixed by the YM coupling parameter [2T], Mo > is an arbitrary constant and A 
the cosmological constant. This infinite mass configuration generalises to higher dimensions the d = 4 
magnetic Reissner-Nordstrom black hole and has a number of interesting properties which are discussed in 
[88j . The generic solutions of the F 2 EYM model have a more complicated pattern (including particle-like 
configuration with a regular origin), but always approach asymptotically the line element (l25|) (see Figure 
!)• 




Figure 2: Typical higher dimensional particle-like (left) and black hole (right) finite mass solutions of the p = 1,2 
EYM theory in asymptotically flat and anti-de Sitter spacetimes. 

The nonexistence result on finite mass solutions is circumvented by adding the appropriate p-YM term(s) 
to the matter Lagrangian. As a result, the EYM system presents (at least) one more coupling constant 
ct 2 = \/ T i / K i T 2i which usually implies a rich structure of the solutions. Various examples were studied: 

• in [18] the particle-like solutions of the system consisting of 1— and 2— Einstein-terms (i.e., the Einstein- 
Gauss-Bonnet system), 1— YM and 2— YM terms in spacetime dimensions d = 6,7,8, thus exhibiting 
a dimensionful constant. Although a pure gauge configuration is approached in the far field, these 
solutions however are not quite direct analogues of the Bartnik-McKinnon solutions because of the 
presence of the dimensionful t 2 constant in the Lagrangian. This is analogous with the the gravitating 
monopole |7J, where a dimensionful constant is also involved. Unlike the latter, however, there were 
no radial excitations in this case. Moreover, the Gauss-Bonnet in the gravity action does not lead to 
any new qualitative features of the solutions, 

• in [20] for particle like solutions of the system consisting of 1— Einstein, 1— YM and 2— YM subsystems 
as above, but in spacetime dimensions d = 5. In addition, in 2D], asymptotically flat black hole 
solutions are constructed. The fixed point properties in d = 5 solutions however differ substantially 
from those of the d — 6, 7, 8 solutions for the same model (see Figure 3). 

symmetries, but modifying the charges in order to take into account the presence of matter fields. A similar approach has been 
used in ref. 21; to assign a finite mass to d > 4 EYM solutions in a F 2 theory with a negative cosmological constant. 
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Figure 3: The value N m of the minimum of the metric function N(r), the mass parameter M as well as the value 
of the metric function a at the origin, er(0), are shown for d = 5, d = 7 asymptotically flat and anti-de Sitter particle 
solutions of the p = 1, 2 EYM theory, as functions of the coupling parameter a 2 — r'f / k\ T2 and several values of 
the cosmological constant A. 

• The fixed point analysis for this model is carried out in [19] . where it is found that in addition to 
the Reissner-Nordstrom fixed point, a new type of fixed point appears. While the Reissner-Nordstrom 
fixed point is typified by the value of the function w(r) = 0, the new type of fixed point is typified by 
the value of the function w(r) = 1 and is referred to as a conical fixed point in [19) . It is further shown 
in [TH], by extending the model judiciously for higher values od d (always keeping only the 1— Einstein 
terms) by higher p YM terms, that this conical singularity appears modulo every 4p dimensions. 

• EYM systems with negative cosmological constant in higher dimensions are also studied in |21| . The 
finite energy solutions in these models exhibit all the properties seen in [18, 20 19 for A = 0. Different 
from the d = 4 case, the higher dimensional AdS solutions necessarily have w(r) — > —1 as r — > oo, i.e. a 
pure gauge configuration is approached in the far field (see Figure 2). As a consequence, asymptotically 
AdS solutions with negative comological constant cannot support nonvanishing A$ solutions, a per the 
argument given in the last but one item in section 2.4. 

• Higher dimensional EYM systems with positive cosmological constant are studied in 22 1. The presence 
of a cosmological horizon leads to a more complicated pattern, where again a conical fixed point appears 
for d = 5. 

All the above listed EYM solutions pertain to models motivated entirely by the criterion of satisfying the 
scaling requirement for finite energy. In a further work [23] in d = 4p, the models were chosen according to the 
criterion that only p— Einstein and p-YM terms appear in the Lagrangiarj^. In these cases no dimensionful 
constant appears in the Lagrangian and the properties of the solutions are entirely similar to those of the 
Bartnik-McKinnon solution (in particular the existence of radial excitations) , except that qualitative features 
are appreciably magnified with increasing p (see Figure 4) . 

There is also the question of non-spherically symmetric EYM solutions in higher dimensions. Since 
all EYM solutions are constructed numerically, the problem here is to relax this symmetry such that the 
numerical process remains tractable. The most s traightforward step would be the imposition of axial 
symmetry in the D— spacelike dimensions, i.e., by imposing spherical symmetry in the (D — 1)— dimensional 
subspace, thus reducing the problem to a 2— dimensional PDE. This can be done readily for arbitrary 
d = D + 1, but unfortunately the implementation of the numerical integration becomes problematic when 
removing the gauge arbitrariness. Instead, for d = 5, a system of 2— dimensional PDE's can be obtained 

8 The ref. 1231 presents also an exact solution for the p-th Einstein- Yang-Mills system in d = 2p + 1 dimensions. 
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Figure 4: The profiles of the metric functions N(r), a(r) and gauge function w(r) are presented for k— node globally 
regular solutions of the p = 2 gravity- Yang-Mills model in d = 8 dimensions. 

when an azimuthal symmetry is imposed in each of the two planes of the 4— dimensional t = const, spacelike 
subspace. In principle, this can be generalized for any odd, d = In + 1, dimensional spacetime, and then 
the reduced problem will be that of a n— dimensional PDE's. This limits one to the bi-azimuthal regime in 
d = 4 + 1, for practical reasons. 

Other than this static result, there are two other indirect d > 4 results in the literature which are not 
spherically symmetric. One is the case where there is a rotation in the two spacelike sub-planes in the 4 + 1 
dimensional case, and the other concerns a rather different topology of the spacetime. 

• Ref. [2S] discussed static solutions in d = 4+1 dimensions of a EYM system with bi-azimuthal symmetry 
in four spacelike dimensions. They generalise the configurations in [18] . [2D] , both particle-like and 
black hole solutions being found to exist. It is interesting that the fixed point structures discovered 
in the spherically symmetric cases in [18l [20l [19j [21j [22] in d = 5, manifest themselves for these bi- 
azimuthally symmetric solutions, although a rigorous fixed point analysis like in [19j is not analytically 
accessible in this case. 

• d = 5 rotating EYM black hole solutions in the usual EYM model (i.e., with a F 2 term only) with 
negative cosmological constant were constructed in 84 . The rotation in question was that of two 
equal angular momenta in the two spacelike 2-planes, which is known to reduce the 2— dimensional 
PDE problem to a 1— dimensional ODE one (i.e. the angular dependence is factorised in the ansatz). 
In this respect, it is a system defined by a single radial variable, but does not not describe a spherically 
symmetric field configuration. Different from the static case, no spinning regular solution is found for 
a vanishing event horizon radius. As expected, the mass of these solutions as defined in the usual 
way, diverges. However, a finite mass can be assigned by using a suitable version of the bounday 
counterterm regularization method 85] • (We would of course expect finite energy solutions in this case 
too, had the F A YM term been included.) 

• In AdS spacetime the topology of the horizon of a black hole solution is no longer restricted to be 
spherical. It can be planar, or, hyperbolic instead. A surprising result reported in [83j is that, for 
A < 0, there are d > 4 asymptotically AdS, finite mass black hole solutions with a planar topology 
of the event horizon, even in a theory without higher derivative terms in the YM curvature. This 
contrasts with the corresponding black holes with a spherical topology of the event horizon, which 
have infinite mass [21] . The case of a hyperbolic topology of the horizon has not been considered yet 
in the literature for d > 4. As in the previous example, the (consistent) Ansatz used for the YM field 
does reduce the PDE's to a system of one dimensional ODE's in terms of a radial variable, but likewise 
does not describe a spherically symmetric field configuration. 
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Figure 5: Fundamental solutions corresponding to several values of the coupling constant a — ti(t p / K p ) 1 ' 2p axe 
shown for p = 4 higher dimensional gravitating YMH monopoles in [24]; the dotted curves are for the monopole in 
flat space. 

In addition to the higher dimensional gravitating YM fields described above, there has been some work 
also studying gravitating monopoles in higher dimensions. Recently a very particular hierarchy of selfgravi- 
tating YMH models in 4p dimensions was studied in detail in [21]. This family of models, whose flat space 
monopoles were constructed in [38], is the most direct generalisation of the d = 3 + 1 Georgi-Glashow model 
(in the BPS limit), and the only one for which the Bogomol'nyi inequalities can be saturated. Both regular 
and black hole solutions have been constructed in [24], which exhibit all the generic properties of the well 
known d — 4 gravitating monopoles (the profiles of typical solutions are exhibited in Figure 5). In higher 
dimensions, there are many other types of monopoles, e.g., that in d = 4 + 1 in [50j . and that in d = 3 + 1 
in [89]. The selfgravitating versions of these are not studied to date. 

Finally, let us mention the case of d > 4 non-Abelian solutions with codimensionfl As mentioned already, 
the situation of d = 5 with one codimension is the only case discussed in a systematic way in the literature, 
mainly for a metric ansatz which is spherically symmetric in a four dimensional perspective 

ds 2 = e -20(r)/v^ + r 2 dn 2^ _ iV ( r ) CT 2 ( r )^ 2 ^ + e ^W/V3 ' ((fa 5)2 j (26) 

x 5 being the extra-direction and <p(r) corresponding to a dilaton field. 

• The KK theory possesses in this case a variety of interesting non-Abelian configurations, including ax- 
ially symmetric generalizations |14j , [15] , [92] , [93] . After performing a KK reduction, they correspond 
to d = 4 particle like and black hole solutions in a Einstein- Yang-Mills-Higgs-U(l)-dilaton theory [94] . 

• EYM black strings and vortices with a cosmological constant were discussed in ref. |95j . 

• The inclusion of of higher order terms of the YM curvature is optional for d — 5 black strings and 
vortices, since they possess a finite mass per unit lenght of the extradimension already in the usual F 2 
theory. In this case, the higher derivative terms do not affect the basic properties of the solutions. 

3.2 EYM solutions with Euclidean signature 

While this review concerns primarily fully gravitating EYM and EYMH solutions in Lorentzian signature, 
it is reasonable to allude to EYM solutions in Euclidean signature, especially since these were the first such 
solutions that appeared in the literature. Such solutions, including those on (Euclidean) Schwarzschild and 

9 A detailed review of these solutions is presented in 1911 . I90| . 
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de Sitter [29, 301 EH] an d Taub-NUT 96[] backgrounds have been studied long ago. More recently there have 
been further investigations [57] [58, 99 of YM fields on fixed backgrounds. 

All of the known EYM fields with Aq ^ and Euclidean signature in the literature are given on fixed 
gravitational backgrounds. None with gravity backreacting on the YM fields is known, which puts these 
solutions on a different footing to those with Lorenzian signature duscussed above. 

Starting with backgrounds which are analytic continuations of relevant solutions with Lorenzian signa- 
tures, we note that 

• The exact solutions of Charap and Duff [29] [30] , as also their higher dimensional counterparts [34] , 
are by construction given on gravitational backgrounds for which the the 2p— form Riemann curvature 
is double-self-dual. Such metrics satisfy the hierarchy of vacuum Einstein equations (with or withour 
cosmological constant) so that by construction, these are EYM solutions on backgrounds of fixed 
curvature. 

• In an effort to go away from fixed backgrounds, a direct numerical method was employed in [100] , using a 
(Euclidean) Schwarzschild metric and static YM fields in d — 4. But being a static field configuration, 
the 'electric' YM potential A assumed the role of a Higgs field and the resulting solutions turned 
out to be self-dual 'deformed Prasad-Sommerfield monopoles', again on a fixed background. These 
solutions are different from those in [29] [30], as shown e.g. by a computation of their action. The 
higher dimensional analogues of this type of solutions in d = Ap were constructed in [34j . 

• In a further development beyond [100J, a number of other static spherically symmetric d — 4 metric 
backgrounds were employed in [101] to construct Euclidean non-Abelian solutions. All these resulted 
in selfdual YM solutions, on the basis of which it was conjectured that for any d = 4 (Euclidean) static 
spherically symmetric metric, the solutions satisfy the d = 4 Yang-Mills self-duality equations. An 
analytic proof of this conjecture has been given in [34] . where the 4p— dimensional analogues of these 
were also constructed, satisfying the self-duality equations l|lip. 

The solutions described in the above three items generalise both the Charap-Duff [3D] and the 
solutions in [101] in d — 4, to d = 4p, the properties of the four dimensional case being generic. 

• Another property of the known gravitating instantons which is related to the fact they are given on 
fixed backgrounds is, that they are always (Euclidean) time independent. Even when an explicit time 
dependence is buit in to the YM Ansatz, it turns out that the solutions are either time indepen- 
den [102 , or in the presence of a cosmological constant, that the Pontryagin charge of the instanton 
is noninteger [103] , We believe that this is a result of having used a static metric. Relaxing this last 
property may be interesting but promises to lead to a nontrivial numerical problem. 

The metric Ansatz ([T7| in all above described Euclidean EYM fields, makes a distinction between the 
(Euclidean) time and the space coordinates. A different type of solutions were found in another setting, the 
metric Ansatz being spherically symmetric in d dimensions [1041 134] 



where f(p) is a function fixed by the gravity-matter field equations, p being the radial coordinate, p = \Z\xn\ 2 
and x^ = x^/p is the unit radius vector. The YM ansatz compatible with the symmetries of the above line 
element is 



where the spin matrices are precisely those used in (fl9l) . ([20l) . 

The resulting reduced one dimensional YM Lagrangian for the p-th term in the YM hierarchy read 



d s 2 = dp 2 + f 2 (p)dnf d ^ 1) , 



(27) 




(28) 
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For any choice of the metric function f(p), the solution of the YM self-duality equation (jTTJ) in d = 4p 
dimensions reads 



w(p) 



1 + c e 



1 - c e 




(30) 



where Co is an arbitrary positive constant. As discussed in [34] , the action of these non-Abelian solutions 
is finite for any value of p. For f{p) — p one recovers the d = 4p generalisation of the BPST instanton 
first found in [36], with w — (p 2 — c)/{p 2 + c). An AdS background f(p) = posinhp/po leads to a d = Ap 
generalisation of the d — 4 AdS selfdual instantons in [104) . with w = (tanh 2 (p/2po) — c)/(tanh 2 (p/2po) + c). 
The d — Ap selfdual instantons on a sphere (euclideanised dS space) are found by taking po — » ipo in the 
corresponding AdS relations. 

As selfdual solutions, these are also fixed background YM fields. Moreover, on curved backgrounds the 
selfduality equation saturating the inequality © can be solved. Thus in Ref. [Mj systems consisting of the 
superposition of two members of the YM hierarchy, say those labeled by p and q, with d = 2(p + q) were 
considered as well. Solutions of such selfduality equations were also discussed in |105j , [106] . 

4 Summary and outlook 

We have reviewed a number of results on Einstein- Yang-Mills (EYM) solutions, with special emphasis on 
the new features that arise for spacetime dimensions d > 4, mainly in Lorentzian signature. The EYM 
solutions that we have considered are those for fully backreacting gravity with matter. It turns out that 
these solutions are all constructed numerically and no relevant closed form solutions are known. As such, 
the question of imposition of symmetries with the aim of reducing the dimensionality of the Euler-Lagrange 
equations becomes a very important feature of these investigations, at least as important as in the case of 
solutions that can be expressed in closed form. 

Higher dimensional EYM fields with Euclidean signature are also mentioned, but only in passing since 
they are not on the same footing as their Minkowskian counterparts. They are exclusively YM fields on fixed 
backgrounds, in all dimensions. 

A salient feature of higher dimensional non-Abelian solutions is that in all d = D + 1 dimensions, for 
D > 4 the usual EYM system cannot support asymptotically flat, finite energy solutions. This is because of 
the inappropriate scaling properties of that system and is remedied by the addition of higher order curvature 
terms. The higher order terms in question are exculsively ones that are constructed with higher order YM 
curvature forms, which in our nomenclature are the p— YM members of the YM hierarchy, the 1— YM being 
the usual YM system. The resulting YM equations contain no higher derivatives of the gauge potential than 
second. The higher order gravitational curvature systems (e.g., Gauss-Bonnet and higher) are not possible 
to exploit for this purpose. Higher order gravities have nonetheless been employed in some contexts, not out 
of necessity, but for emphasising qualitative features of certain EYM and EYM-Higgs (EYMH) solutions, 
which get magnified if dimensions of the consituent terms in the system are suitably matched (see e.g. [23 ). 

The necessity of employing higher curvature members of the YM hierarchy to enable finite energy holds 
both for asymptoticall flat EYM, as well as in the presence of a cosmological constant. There is however an 
exception to this rule, namely in the case of asymptotically AdS EYM black hole solutions with a Ricci flat 
horizon geometry. In that case, which is of particular interest for applications to AdS/CFT, it turns out that 
the usual EYM system can support finite energy solutions in all dimensions [83]. There, the nonexistence 
proof for solutions with nonvanishing 'electric' YM potential does not hold. 

Concerning the physical justification of employing higher order YM and Riemann curvature terms, one 
notes that these occur in the low energy effective action of string theory [TBI [17]. Indeed there is some 
controversy on the precise structure of such terms, especially in the YM case [107111081 1109] . but we do not 
take account of these considerations here. From the point of view of applications for higher dimensional 
EYM, these can be used to extend the results of [110111111,1112] to higher dimensions. Some of these authors, 
and |113[ 1114] , employ only solutions in closed form, but further developments would necessitate the fully 
gravitating solutions which are evaluated only numerically. In any case, our focus here was exclusively on 
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the existence and the generic properties of EYM solutions in higher dimensions, rather than their physical 
applications. 

Nearly all gravitating non-Abelian solutions in d = D + 1 dimensions with D > 4 reviewed in this work 
are static and spherically symmetric. This contrasts with the situation in 3 + 1 spacetime dimensions where 
axial symmetry in 3 space dimensions is really azimuthal symmetry and the non-Abelian field configurations 
are encoded with a vortex (winding) number. The latter turns out to be an essential tool in the numerical 
constructions |251 126] . In higher than four spacetime dimensions however, axial symmetry implies spherical 
symmetry in one dimension lower than D and there is no winding number associated with the axially 
symmetric fields. However, the removal of the gauge arbitrariness turns out to be a much harder problem in 
this case, presenting a technical obstacle. Thus, the only non-spherically symmetric EYM solutions studied 
to date are ones in 4+ 1 spacetime with bi-azimuthal symmetry, when there are two vortex numbers encoding 
the symmetries of the field configurations [28] . 

Looking further ahead we should note that in recent years it has become clear that as the dimension d 
increases, the phase structure of the (non-spherically symmetric) solutions of the Einstein equations becomes 
increasingly intricate and diverse, already in the vacuum case (see e.g. the recent work |115j ). It is very likely 
that, given the interplay in this case between internal group symmetries and the spacetime symmetries, the 
extension of known such vacuum solutions to a non-Abelian matter content would lead to a variety of new 
unexpected configurations. This is a promising but technically difficult direction for the future. 

Finally we mention the higher dimensional Euclidean EYM fields, which is of marginal interest here since 
none of the known such solutions are genuinely selfgravitating, but rather YM fields on fixed backgrounds. 
In all dimensions, including four, these turn out to be selfdual YM fields, and as such are restricted to 
even dimensions only. Another aspect of this restriction turns out to be that these YM fields appear to 
be (Euclidean-) time independent and are not YM instantons at all. It would be interesting to construct 
non-selfdual Euclidean EYM fields and inquire whether these, if they exist, describe genuine time dependent 
instantons. 
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